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In this paper we present a detailed analysis of the structural, electronic, and optical properties 
of an m-plane (In,Ga)N/GaN quantum well structure grown by metal organic vapor phase epitaxy. 

The sample has been structurally characterized by x-ray diffraction, scanning transmission electron 
microscopy, and 3D atom probe tomography. The optical properties of the sample have been 
studied by photoluminescence (PL), time-resolved PL spectroscopy, and polarized PL excitation 
spectroscopy. The PL spectrum consisted of a very broad PL line with a high degree of optical 
linear polarization. To understand the optical properties we have performed atomistic tight-binding 
calculations, and based on our initial atom probe tomography data, the model includes the effects of 
strain and built-in field variations arising from random alloy fluctuations. Furthermore, we included 
Goulomb effects in the calculations. Our microscopic theoretical description reveals strong hole 
wave function localization effects due to random alloy fluctuations, resulting in strong variations in 
ground state energies and consequently the corresponding transition energies. This is consistent with 
the experimentally observed broad PL peak. Furthermore, when including Goulomb contributions 
in the calculations we And strong exciton localization effects which explain the form of the PL 
decay transients. Additionally, the theoretical results confirm the experimentally observed high 
degree of optical linear polarization. Overall, the theoretical data are in very good agreement with 
the experimental findings, highlighting the strong impact of the microscopic alloy structure on the 
optoelectronic properties of these systems. 

I. INTRODUCTION 


(In,Ga)N/GaN quantum well (QW) structures grown 
on polar c-plane substrates are the building blocks for 
nitride-based light emitters operating in the blue spec¬ 
tral regioniii^ However, when trying to push the emission 
wavelength into the green to yellow spectral region the 
device performance is hampered, among other factors, 
by the presence of strong electrostatic built-in fields, aris¬ 
ing in part from the strain-dependent piezoelectric polar¬ 
ization and in part from the strain-independent sponta¬ 
neous polarizationi^ One consequence of the built-in held 
is the so-called quantum conflned Stark effect (QGSE) 
which results in a shift of the emission to longer wave- 
lengthi^i^ Also the intrinsic held leads to a spatial sep¬ 
aration of electron and hole wave functions, causing a 
reduced radiative recombination rate;^ an effect that 
can be particularly undesirable for high-efhciency opto¬ 


electronic devices. To circumvent these effects arising 
from the intrinsic built-in helds, which fundamentally 
are caused by the growth along the polar c axis, signif¬ 
icant research has been directed towards the fabrication 
of semi- and nonpolar structuresIn the case of semi- 
and nonpolar planes the c axis is at a nonvanishing an¬ 
gle with respect to the growth direction. In semipolar 
planes residual built-in helds are expected and observed, 
even though these helds are signihcantly reduced In 
terms of built-in held reduction, nonpolar QW systems 
are ideal since in these structures the c-axis lies in the 
growth plane. Thus, in a perfect nonpolar system, there 
is no discontinuity in the built-in polarization vector held 
at the heterointerfaces, leading ideally to a held-free sys- 
temiii^ Consequently, in such a structure the radiative 
recombination rate should be much higher than in polar 
QW structures. Additionally, (In,Ga)N QWs grown on 
nonpolar planes offer the possibility of acting as highly 
efficient sources of linearly polarized light /^d^d^ which 
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may be of practical use in, e.g., back-lit liquid crystal dis¬ 
plays.— The potentially high degree of optical linear po¬ 
larization (DOLP) in a nonpolar QW originates mainly 
from differences in the effective hole masses along the 
growth direction, which leads to a lifting of the degener¬ 
acy of the highest lying p-like valence bands i^^i^^ 

However, the detailed consequences of random al¬ 
loy fluctuations on the ideally zero built-in held and 
the optoelectronic properties of nonpolar systems have 
been largely ignored. Experimental studies on c-plane 
(In,Ga)N/GaN QW systems reveal that random alloy 
huctuations lead to carrier localization effects, which can 
dominate the electronic and optical properties of these 
structures— *21 Clear experimental indications of the im¬ 
portance of alloy fluctuations on the properties of non¬ 
polar (In,Ga)N QWs have also been presented.— To 
fully understand the structural, electronic, and optical 
properties of nonpolar (In,Ga)N QWs requires the use of 
a range of advanced experimental and theoretical tech¬ 
niques. For example, 3D atom probe tomography (APT) 
can be used to help to identify whether the QW should 
be described by a clustered or random alloy— This 
structural information can then form the basis of de¬ 
tailed theoretical modeling. However, from a theoretical 
viewpoint, polar and nonpolar (In,Ga)N QWs are widely 
treated as homogeneous systems described by average 
material parameters. Obviously, these approaches com¬ 
pletely neglect wave function localization arising from ef¬ 
fects that could be attributed to the microscopic alloy 
structure. Recently, continuum-based approaches have 
been modihed to include random alloy fluctuations.—^— 
Even though such an approach captures some of the lo¬ 
calization effects introduced by alloy fluctuations it can¬ 
not reveal the microscopic origin of localization effects, 
including in particular the presence of In-N-In-N chains 
as shown by density functional theory (DFT)— Un¬ 
fortunately, realistic QW systems cannot be treated in 
the framework of DFT due to the very large number of 
atoms (> 10^) required in the calculations. Additionally, 
especially for nonpolar systems where the strong macro¬ 
scopic electrostatic built-in field is absent, DFT calcula¬ 
tions would have to be coupled with solving the Bethe- 
Salpeter equation or time-dependent DFT would have 
to be used to account for excitonic effects. These effects 
could be very important when comparing the results with 
photoluminescence (PL) measurements.— To be able to 
take all these factors into account, namely the large num¬ 
ber of atoms, the atomistic description of the electronic 
structure, plus Goulomb correlation effects, semiempiri- 
cal approaches, such as tight-binding^^— (TB) or em¬ 
pirical pseudopotential methods^i^ coupled with con¬ 
figuration interaction (Cl) schemes^i^li^ are required. 

In this work we combine advanced experimental 
and theoretical methods to analyze the structural, 
electronic and optical properties of nonpolar m-plane 
(In,Ga)N/GaN QWs. Our approach allows us to clar¬ 
ify the impact of alloy fluctuations on the properties 
of nonpolar (In,Ga)N structures. We have conducted a 


comparative study of the predicted and measured optical 
properties of a 2 nm thick m-plane (In,Ga)N/GaN QW 
structure grown on a freestanding GaN substrate. The 
heterostructure was grown by metal organic vapor phase 
epitaxy (MOVPE) and characterized by x-ray diffrac¬ 
tion (XRD), scanning transmission electron microscopy 
(STEM), and 3D APT measurements. The detailed opti¬ 
cal characterization of the sample was carried out by PL, 
polarized PL excitation (P-PLE) spectroscopy, and PL 
time decay measurements. The measured optical prop¬ 
erties are compared with theoretical results. The the¬ 
oretical framework is based on an atomistic TB model 
coupled with GI calculations to account for Goulomb ef¬ 
fects. Our microscopic model takes input from experi¬ 
mentally determined structural properties and accounts 
for local strain and built-in field fluctuations arising from 
alloy fluctuations in the QW region. 

The TB model assumes that the In atom distribu¬ 
tion inside the QW is close to a random alloy, and we 
briefly discuss this assumption in the context of available 
APT data. The optical characterization reveals strong 
Stokes shifts and broad PL linewidth, which is indicative 
of strong carrier localization effects. Additionally, and in 
stark contrast to c-plane systems, our time-resolved mea¬ 
surements exhibit decay transients that are single expo¬ 
nential. Also, we hnd in general an extremely high DOLP 
(> 90%). Using the experimentally obtained structural 
information as input for atomistic TB calculations we 
find very good agreement between theory and the opti¬ 
cal spectroscopy measurements. 

The theoretical analysis reveals strong hole wave func¬ 
tion localization, which leads to a broad emission spec¬ 
trum. When we include the Goulomb interaction be¬ 
tween the electrons and holes in our calculations we find 
that electron wave functions localize about the hole wave 
functions, giving rise to localized excitons. This find¬ 
ing is consistent with the experimentally observed form 
of the decay transients. Also, the calculated PL transi¬ 
tion energy is in good agreement with the experimentally 
measured value. Additionally, analyzing the wave func¬ 
tion character of the valence band edge (VBE) shows a 
high DOLP similar to the experimental data. Overall, 
this combined experimental and theoretical analysis pro¬ 
vides clear insights into the basic physical properties of 
TO-plane (In,Ga)N/GaN QWs. 

The paper is organized as follows. The details of the 
growth of the sample are given in Sec.|TTl This is followed 
by the description of the applied structural characteriza¬ 
tion techniques, the results obtained, and how they are 
implemented in the theoretical framework (Sec. IBII) . In 
Secs. lIV Al and lTV^ we describe the setup and the results 
of the experimental optical characterization of the m- 
plane (In,Ga)N/GaN QWs, respectively. In Sec. HV CI the 
ingredients of our theoretical framework are discussed, 
while Sec. IIVDI presents the theoretical results and the 
detailed comparison with the experimental data. Finally, 
we summarize our results in Sec. 
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II. SAMPLE GROWTH DETAILS 

The TO-plane (In,Ga)N/GaN multiple-QW structure 
studied was grown by MOVPE in a Thomas Swan 6x2 
in close-coupled showerhead reactor on a freestanding m- 
plane (1-100) GaN substrate. The misorientation of the 
substrate is 2° ±0.2° in the negative c direction. The 
substrate, provided by Ammono S. A.,— has a negligible 
basal plane stacking fault density and a threading dislo¬ 
cation density of < 5x10"^ cm“^. Onto this freestanding 
substrate an initially undoped 2 thick homoepitax- 
ial GaN buffer layer was grown. For this growth step 
trimethylgallium (TMG) and ammonia (NH 3 ), with H 2 
as the carrier gas, were used. A reactor pressure of 100 
Torr at 1050°G was applied. Subsequently five (In,Ga)N 
QWs were grown with a quasi-2T method, discussed in 
detail elsewhere.— For the QW growth a reactor pressure 
of 300 Torr was applied. TMG, trimethylindiuin (TMI), 
and NH 3 have been used as precursors, with N 2 as the 
carrier gas. The (In,Ga)N QWs were grown at a temper¬ 
ature of 735°G. After the growth of each QW period a 1 
nm GaN cap was grown at the (In,Ga)N growth temper¬ 
ature. The flows of TMG and NH 3 were kept constant 
during the 90 s temperature ramp to 860°G. After this 
90 s temperature ramp the growth of the barrier material 
was completed. 


III. STRUCTURAL CHARACTERIZATION: 

SETUP, RESULTS AND INPUT INTO 
THEORETICAL FRAMEWORK 

Having presented growth details of the sample, we 
now turn to the analysis of its structural properties. In 
Sec. IIII Al we briefly introduce the experimental tech¬ 
niques applied to analyze the structural properties of the 
TO-plane sample and the outcome of these studies. We 
discuss in Sec. ME how the experimental results serve 
as input for theoretical supercell (SG) structures, which 
underlie the atomistic calculations of the electronic and 
optical properties of the sample in question. 


A. Structural characterization: Experimental setup 
and results 

To analyze the structural properties of the TO-plane 
(In,Ga)N/GaN QW system, different techniques have 
been applied to obtain a comprehensive picture. The 
Ina;Gai_a;N QW and barrier widths as well as the In 
fraction (cc) of the sample were determined by XRD. The 
sample was analyzed using a high-resolution x-ray diffrac¬ 
tion MRD diffractometer from Panalytical equipped with 
a symmetric 4-bounce monochromator and a triple-axis 
analyzer to select the GuKal wavelength, oj — 29 scans 
of the brightest symmetric reflection ( 1 - 100 ) were used 
to analyze the multiple QWs, with a range of 10 de¬ 
grees in w. The measured QW width was 2.0±0.3 nm 



FIG. 1: APT analysis of the first QW (closest to the sub¬ 
strate); box size is 55 x 55 x 11 nm®. (a) 3D reconstruction 
of the first QW showing 30% Ga atoms and 50% In atoms; 
(b) upper and (c) lower interfaces of the QW rendered as 
isosurfaces at 5% In content. 


with a barrier width of 6.1±0.3 nm. The In composition 
was determined to be 18.3±0.7%. The sample was also 
characterized with an FEI Osiris fitted with an extreme- 
Schottky field emission gun operated at 200 keV. The 
measured QW and barrier widths confirmed the XRD 
measurements. 

APT experiments were conducted in pulsed laser mode 
with a pulse energy of 0.012 nJ using a Cameca LEAP 
3000X HR, where the base temperature of the sample 
was set at ~30 K and a constant detection rate of 0.01 
atoms per pulse was employed. The IVAS software pack¬ 
age (GAMECA Version 3.6.6) was used to carry out APT 
reconstruction and analysis was informed by the thick¬ 
nesses of the (In,Ga)N and GaN layers measured by XRD 
and transmission electron microscopy (TEM) and also 
the geometry of the tip as measured by scanning elec¬ 
tron microscopy (SEM). The In content of the sample 
was analyzed using concentration profiles through each 
QW computed using a proximity histogram (proxigrain) 
approach as described in Ref. l26l . The mean In fractions 
across the five QWs based on the maxima of the In con¬ 
centration profiles is 17.0±0.6%, which is in good agree¬ 
ment with the XRD data. While frequency distribution 
analysis of c-plane QWs provides no evidence for the for¬ 
mation of nonrandom indium clusters;^ recent studies of 
nonpolar a-plane QWs with a similar In content (15%) 
to the TO-plane samples studied here have presented clear 
evidence for the presence of clustering— For the sample 
under discussion here, the question of clustering is less 
clear-cut. Evidence of clustering is sensitively dependent 
on the exact parameters of the analysis used, and the is¬ 
sue is still under investigation. However, in this paper we 
model the QW as a random alloy and assess the extent to 
which this model can accurately reproduce the observed 
optical properties, which can provide an insight into the 
relative importance of slight deviations from randomness 
in these QWs, if present. More details of the APT anal- 
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ysis will be reported in Ref. |4l|. 

An isoconcentration surface analysis of the APT data 
similar to that performed by Galtrey et al^ has been 
performed to assess whether the roughness of the up¬ 
per interface of the QWs [at which GaN is grown on 
(In,Ga)N] is rougher than the lower interface of the QWs 
[at which (In,Ga)N is grown on GaN] and typical data are 
shown in Fig.[T]for an indium concentration at the isosur¬ 
face of 5%. This yields a roughness for the lower surface 
of 0.24 nm (root mean square), but a larger roughness for 
the upper surface of 0.48 nm, with the increased rough¬ 
ness largely relating to raised island-like features a few 
nm in diameter and 1-2 monolayers in height, consis¬ 
tent with the monolayer and bilayer well-width fluctua¬ 
tions previously observed on the c plane. High-resolution 
scanning transmission electron microscopy imaging of the 
surfaces of the QWs also indicates the existence of sur¬ 
face roughness. Although the projection problem makes 
it difficult to distinguish the exact nature of the rough¬ 
ness observed in STEM, the data are consistent with the 
APT observations. 


B. Structural characterization: Input into the 
theoretical framework 

Based on the experimental data discussed in Sec. lIII Al 
we take the m-plane QW width to be approximately 2 
nm. We include in our calculations disklike 2 monolayer 
thick well-width fluctuations at one of the QW interfaces. 
This is a reasonable representation of the well-width fluc¬ 
tuations that have been observed in APT and is con¬ 
sistent with earlier approaches to the modeling of such 
featuresj^ However, for the m-plane QW structure con¬ 
sidered here, we do not expect that 2 monolayer thick 
well-width fluctuations are of central importance for the 
description of the electronic and optical properties. This 
assumption is based on the fact that in an m-plane QW 
the macroscopic field is absent. Therefore, in contrast to 
a c-plane system, the wave functions in an m-plane sys¬ 
tem are not localized near the interface between QW and 
barrier material, since there is no field causing this form 
of localization. Gonsequently, the wave functions are ex¬ 
pected to be far less sensitive to the shape and diameter 
of the well-width fluctuations. As we will show below, 
this assumption is consistent with the result of our cal¬ 
culations. The In content in the sample has been mea¬ 
sured to be 18.3±0.7% (XRD) and 17.0±0.6% (APT), 
respectively. For our calculations we have set the QW 
In content to 17%, close to the measured values. We use 
here SCs with a size of 10x9x10 nm^, thus containing ap¬ 
proximately 82,000 atoms, with periodic boundary con¬ 
ditions. Based on the APT data discussed in Sec. IHI Al 
we treat (In,Ga)N as a random alloy. To be able to per¬ 
form a detailed comparison with experiment and to allow 
for reasonable statistical averages, our calculations have 
been repeated 75 times for the fixed In fraction of 17%. 
This allows us to realize a large number of microscopi¬ 


cally different random configurations in a large SG, where 
we have nominally the same In content in the SC but the 
positions of the In atoms have been changed randomly. 

IV. ELECTRONIC AND OPTICAL 

PROPERTIES: EXPERIMENT AND THEORY 

Having investigated the structural properties of the 
sample, we turn now to study the electronic and optical 
properties of the m-plane (In,Ga)N/GaN QW system. 
We start in Sec. IIV Al by describing the experimental 
methods used for the optical characterization, followed 
by the discussion of the experimental results in Sec. IIV Bl 
We then briefly introduce the theoretical framework in 
Sec. IIV Cl while Sec. IIV PI focuses on the comparison of 
theoretical and experimental data. 

A. Experimental setup 

The PL and P-PLE studies were carried out either us¬ 
ing excitation from a cw He/Cd laser or using a combi¬ 
nation of a 300 W xenon lamp and monochromator as a 
fixed or tunable wavelength excitation source. The sam¬ 
ple was mounted in the cryostat so that the c axis of 
the GaN was horizontal. The PL from the sample was 
analyzed by a 0.85 m double-grating spectrometer and a 
Peltier-cooled GaAs photomultiplier using standard lock- 
in detection techniques. The spectral dependence of the 
DOLP of the emission was determined by measuring the 
PL spectra polarized in the plane parallel (/||) and per¬ 
pendicular {I±) to the c axis of the sample. Once the 
spectra are obtained, the DOLP is calculated iiom^ 

DOLP = , (1) 

+ 7 || 

where I± (/||) is the PL intensity for the electric field E 
perpendicular E T c (parallel E || c) to the c axis. 

A time-correlated single-photon counting system was 
used for the time decay measurements. For these stud¬ 
ies the exciting light was generated by a mode-locked 
Tiisapphire laser with a frequency-tripled emission out¬ 
put of 4.88 eV. A spectrometer followed by a micro chan¬ 
nel plate detector was used to detect the emitted light. 

B. Experimental results: Optical characterization 

In this section we describe the results of our PL, P- 
PLE, and time decay measurements at 10 K on the m- 
plane (In,Ga)N/GaN QW sample described above. Fig¬ 
ure [2] shows the P-PLE spectra with the plane of polar¬ 
ization of the incident light either perpendicular (E T c) 
or parallel (E || c) to the c axis of the crystal. Figure [D 
also displays the (unpolarized) PL spectrum of the sam¬ 
ple (solid black line). We note several specific aspects 
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FIG. 2: Low-temperature (10 K) PL and PLE spectra for the 
m-plane (In,Ga)N/GaN QW structure under consideration. 
The PLE spectra are shown with the plane of polarization of 
the incident light perpendicular (blue dashed) or parallel (red 
dashed-dotted) to the c axis of the crystal. The PL spectrum 
is given by the black solid line. For PL and PLE measure¬ 
ments the excitation spot was 0.02 cm^ with an excitation 
power density of ~ 10 /rW/cm^. 



FIG. 3: (Color online) The unpolarized PL spectrum is dis¬ 
played by the solid black line. The data are taken under 
time-averaged conditions, where the average power density of 
the laser source was 1 W/cm^ using 100 fs pulses at a repeti¬ 
tion rate of 800 kHz and a spot size of 2 x 10~® cm^. The inset 
shows the time-dependent PL results. The detection energy 
for the PL transients was 2.934 eV. The spectral dependence 
of corresponding time constants r is given (blue squares). All 
measurements are performed at low temperatures (10 K). 


of the PL spectrum. First, the full width at half max¬ 
imum (FWHM) of the PL spectrum is 135 meV. This 
value is much larger than typical FWHM values found 
in c-plane systems4^ Although it should be noted that 
the PL spectrum exhibits an extended low-energy tail, 
we have shown previouslj^ that in our samples this low- 
energy tail arises from recombination in QWs present on 
semipolar facets which form at step bunches associated 
with the miscut of the GaN substrate. These semipo¬ 
lar QWs have higher QW width and In content than the 
nonpolar QWs on the adjoining m-plane terraces on ei¬ 
ther side of the step, and thus emit at longer wavelength. 
We do not aim to include such structures in our model 
or reproduce this aspect of the spectra here. The emis¬ 
sion from the semi-polar QWs is of secondary importance 
in determining the magnitude of the spectral width. Sec¬ 
ond, the energy difference (Stokes shift) between the peak 
of the lowest exciton transition in the P-PLE spectrum 
and the peak of the PL spectra is 180 meV. We also note 
an energy splitting between the two lowest exciton tran¬ 
sitions for the two polarizations of the excitation light of 
35 meV. We attribute this splitting mainly to the crystal 
field splitting energy in the (In,Ga)N system4^ 

The values of the linewidth and Stokes shift, based on 
the assumption that the emission is an intrinsic process, 
are strongly suggestive of emission involving strongly lo¬ 
calized carriers. As to the nature of the localization we 
anticipate that, similarly to c-plane systems, alloy fluc¬ 
tuations will play a major role. 

To help further understand the nature of the recombi¬ 
nation we performed PL time decay measurements across 
the spectrum. The results of these measurements are de¬ 


picted in Fig. [3l Interestingly, our measurements reveal 
that the PL decays are single exponentials with charac¬ 
teristic time constants r = 350 ± 20 ps over the majority 
of the spectrum. Throughout this work we assume that 
the recombination is purely radiative and the values of 
T are the radiative decay constants. This assumption is 
supported by temperature-dependent measurements (not 
shown) and the fact that the values of r found here are 
similar to those of Marcinkevicius et ah^, also measured 
at low temperature. The data reported by Marcinkevi¬ 
cius et at low temperature were shown to be not 
influenced by nonradiative recombination. Only on the 
low-energy side of the spectra did we find more slowly 
decaying emission as we have reported elsewherei^ The 
emission on the low-energy side is attributed to semipo¬ 
lar QWs at step bunches described above. The decay 
time from the emission associated with these semi-polar 
QWs occurs on a longer time scale due to the locally 
present built-in field, which separates the electron/hole 
wave functions. 

It is important to note that the radiative decay dy¬ 
namics observed here are in stark contrast to c-plane 
(In,Ga)N QW structuresFirst, in polar (In,Ga)N 
QWs the decay occurs over a much longer time scale 
due to the polarization field perpendicular to the plane 
of the QWs4^i^^— Second, the decay curves are non¬ 
exponential due to the variable in-plane separation of 
the separately localized electrons and holes4Si^ Single¬ 
exponential decays in nonpolar QWs have also been re¬ 
ported by other groupsi^i^ The explanation proposed 
for the nature of the decay transients is that recombina¬ 
tion involves localized excitons4^i^ 










6 



Ph 

hJ 

O 

o 


FIG. 4: (Color online) The (unpolarized) PL emission spec¬ 
trum for the m-plane (In,Ga)N/GaN QW sample is shown 
by the black solid line while the spectral dependence of the 
DOLP is given by the red dashed line. All measurements have 
been performed at low temperatures (10 K). The cw excita¬ 
tion power density was 12 W/cm^ with an excitation spot of 
5 X 10“® cm^. 


As discussed in the introduction, nonpolar QWs should 
exhibit a high DOLP. This should be particularly the 
case at 10 K, bearing in mind the splitting of the low¬ 
est valence band states as revealed by the PLE measure¬ 
ments shown in Fig. [2] The measured DOLP of the m- 
plane (In,Ga)N/GaN QW sample at 10 K, based on the 
method described in Sec. IIV Al is shown in Fig. |4] to¬ 
gether with the (unpolarized) PL spectrum. We find a 
very high DOLP (> 90%), which is in good agreement 
with other experimental data on m-plane (In,Ga)N/GaN 
Theoretically, based on continuum-like de¬ 
scriptions, at the VBE one would expect a very high 
DOLP for several reasons. First, assuming growth along 
the y direction, the QW confinement lifts the degener¬ 
acy of the |X)- and |y)-like valence band states due to 
their different effective masses along the growth direc- 
tionj ^^d^i^^ Second, jAT)- and |Z)-like states are energet¬ 
ically separated by the crystal field splitting energy in 
(In,Ga)N alloys4^ Thus, the VBE is expected to be pre¬ 
dominately |Ar)-like. Figure |4] shows that the DOLP is 
almost constant over the entire spectrum. Only on the 
low-energy side of the spectrum do we find a small re¬ 
duction in the DOLP. This reduction is to be expected as 
the low-energy part of the PL spectrum originates from 
recombination involving semipolar QW states. The de¬ 
tailed strain state of such QWs on nanoscale facets asso¬ 
ciated with step bunches is not straightforward to predict 
and might cause increased band mixing effects and con¬ 
sequently a reduced DOLP. A detailed investigation of 
these semipolar contributions is beyond the scope of the 
present study. Overall, given that the PL linewidth ex¬ 
ceeds the valence band splitting, it is not immediately 
obvious why the DOLP, shown in Fig. 01 is virtually con¬ 
stant across the spectrum. One possible answer is that 
the PL emission spectrum originates almost entirely from 


transitions involving localized states in different potential 
fluctuations. Therefore, calculations including random 
alloy fluctuations are required to shed more light on this 
behavior. In the next section we introduce the theoretical 
framework used to address these questions. 


C. Theoretical Framework 

We use an atomistic TB model to study the electronic 
and optical properties of m-plane (In,Ga)N-based QWs. 
This approach allows us to include (random) alloy fluc¬ 
tuations in the QW region and the resulting local strain 
and built-in potential fluctuations on a microscopic level. 
The strain field calculations are based on a modified va¬ 
lence force field model accounting for electrostatic correc¬ 
tions in addition to bond bending, bond stretching and 
related cross-termsi^ Local built-in field fluctuations are 
treated on the basis of our recently developed local polar¬ 
ization theoryi^ Strain and built-in held huctuations are 
included in the sp^ TB model presented in Ref. The 
bulk TB parameters were determined from htting TB 
band structures to hybrid functional (HSE) DFT band 
structures, showing a very good agreement between the 
TB and HSE-DFT band structures. The TB parameters 
and the TB model are discussed in detail in Ref. The 
application of the TB model to QW structures, includ¬ 
ing the coupling of the TB model to Cl calculations to 
include excitonic effects, is discussed in detail in Ref. HI. 
The model has already been benchmarked against experi¬ 
mental and DFT data on bulk (In,Ga)N alloys, revealing 
a very good agreement between experiment, DFT, and 
TB results Moreover, our recent calculations on elec¬ 
tronic and optical properties of c-plane (In,Ga)N/GaN 
QWs show also a very good agreement with available ex¬ 
perimental data<^ The model can directly be applied, 
without modihcation of the theoretical framework, to m- 
plane structures. Here, only the orientation of the QW 
structure in the simulation SC has to be changed com¬ 
pared to a c-plane system. The numerical and structural 
details of the system considered here (QW width, In con¬ 
tent, number of atoms in the SC, SC size, etc.) have 
already been discussed in Sec. IIIIBI 


D. Theoretical results: Electronic structure and 
optical properties 

Equipped with the knowledge about the experimental 
data, we present here the results of the TB calculations. 
Figure[5]shows the variation of the single-particle electron 
AFqs and hole A£'q 3 ground state energies with respect 
to the average ground state energy for electrons Eq^ and 
holes Eq^ , respectively, as a function of the configuration 

(Config) number n ^AAgg = (Conhg n) — 

From Fig. [5] we can infer that variations in the electron 
ground state energies are much smaller than the varia- 
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FIG. 5; (Color online) Variation of the electron AEqq (blue 
square) and hole AEqs (red circle) ground state energies 
with respect to their average ground state energies E^g 
(^AEqq = (Config n) — E^g^j . The results are shown as 
a function of the configuration (Config) number n. 


tions in the hole ground state energies (Ai^^g ^ AEq^) . 
This is also reflected in the calculated standard devia¬ 
tions a. For electrons we And a standard deviation of 
CT® = 9.8 meV, while the standard deviation for the hole 
ground state energy is cr^ = 33.7 meV. We have observed 
a similar behavior in c-plane (In,Ga)N/GaN QWsi^ 

The strong variations in the hole ground state energies 
in c-plane systems are attributed to hole wave function lo¬ 
calization effects due to random alloy fluctuations 
To confirm this behavior in a nonpolar system, Fig. [S] 
shows isosurfaces of the single-particle electron (red) and 
hole (green) ground state charge densities for configu¬ 
rations 5, 20, and 75, which reflect typical situations 
observed in the random alloy as we will discuss below. 
The light and dark isosurfaces correspond to 5% and 
25% of the maximum charge density values, respectively. 
From Fig. [S] one can infer that the random alloy fluc¬ 
tuations lead to strong hole wave function localization 
effects, independently of the configuration. In general we 
find that the alloy fluctuations have a much weaker effect 
on the electron wave functions, when compared with the 
holes. As expected from our discussion in Sec. IIIIBl 2 
monolayer thick well-width fluctuations do not contribute 
significantly to wave function localization effects. Un¬ 
less very deep well-width fluctuations are present, which 
would generate large c-oriented surface areas, we do not 
expect that well-width fluctuations change the localiza¬ 
tion characteristics of electron and hole wave functions 
in comparison to the here presented data. Therefore, 
we conclude that in m-plane QW structures, grown on 
freestanding GaN substrates, localization effects can be 
mainly attributed to alloy fluctuations. Looking at the 
displayed configurations in more detail, for configuration 




FIG. 6: (Color online) Single-particle ground state electron 
(red) and hole (green) charge densities for configurations 5 
(Config 5), 20 (Config 20), and 75 (Config 75). Results are 
shown for different view points. Light (dark) isosurfaces cor¬ 
respond to 5% (25%) of the maximum charge density value. 
Dashed lines indicate the QW interfaces. 



































































FIG. 7: (Color online) Single-particle (black circles) and ex- 
citonic (red squares) transition energies as a function of the 
configuration number n. The average single-particle transi¬ 
tion energies are indicated by the (black) dashed line, while 
the (red) dashed-dotted line indicates the average excitonic 
transition energy. 


5, the electron wave function is almost spread out over 
the entire QW area, with perturbations introduced by the 
alloy fluctuations. We observe slightly stronger perturba¬ 
tions in the case of configurations 20 and 75 (cf. Fig. [6]). 
Overall there is a clear difference in the wave function 
localization characteristics between electrons and holes. 
All this shows that a continuum-based description might 
give, to a first approximation, a reasonable description 
of the ground state electron wave functions but would 
fail for the description of the hole ground states. One 
should also note that the differences in the localization 
features of the different electron single-particle states, as 
discussed above, are important. For instance, in configu¬ 
ration 75 the electron charge density is mainly localized 
in the region where the hole ground state is localized. 
Consequently, one could expect a higher wave function 
overlap between electron and hole. For configuration 5, 
the electron ground state charge density has stronger con¬ 
tributions (dark red isosurfaces) in areas spatially sepa¬ 
rated from the strongly localized hole ground state. This 
is even more pronounced for configuration 20 (cf. Fig. [6]). 

In general, the observed asymmetry in the localiza¬ 
tion characteristics of electron and hole ground state 
wave functions is consistent with DFT-based calculations 
for bulk (In,Ga)N alloysr^*^ and explains the observed 
asymmetry in the electron and hole ground state energy 
variations displayed in Fig. [5l While the local atomic 
arrangement (In-N-In chains) is of secondary importance 
for the electrons, it is of central importance for the hole 
ground states, as shown by DFT-based calculations.—!^ 
Thus, different microscopic arrangements of In atoms 
will give very different hole ground state energies. This 
picture of localized hole and delocalized electrons is not 



FIG. 8: Excitonic binding energy Ex,b as a function of the 
configuration number n. The average excitonic binding en¬ 
ergy Ex,b is indicated by the dashed line. 


quite compatible with the experimentally observed decay 
times being constant across the spectrum, as we would 
still expect some variation in the electron and hole in¬ 
plane wave function overlap leading to a range of ra¬ 
diative recombination times, i.e., nonexponential decay 
curves. It should however be noted that the effect would 
be much less marked than in the case of polar QWs where 
the electrons and holes are localized separately. We will 
discuss this discrepancy later. 

The variations in the single-particle energies will also 
affect the single-particle ground state emission energies. 
This is shown in Fig.[7l The single-particle transition en¬ 
ergies (black circles) show large variations between 2.81 
eV and 3.07 eV. The average single-particle transition 
energy is around 2.93 eV (black dashed line). However, 
to be able to compare the theoretical ground state emis¬ 
sion energies with the experimental PL emission energies, 
we have to include Coulomb (excitonic) effects. We use 
the Cl scheme described in detail in our previous work, 
to include excitonic effects in the description— We as¬ 
sume here a single electron-hole pair. Thus, electron- 
electron and hole-hole Coulomb interactions are not con¬ 
sidered. We neglect electron-hole exchange contributions 
since these are small corrections on the energy scale rel¬ 
evant for the discussion of our results. To describe the 
excitonic many-body wave function we include 5 electron 
and 15 hole states in the Cl expansion. The calculated 
excitonic ground state energies (red squares) are shown 
in Fig. [7] again as a function of the configuration num¬ 
ber n. The average excitonic transition energy is given 
by the (red) dashed dotted line. From Fig. [7] we con¬ 
clude that Coulomb effects introduce a significant shift 
and broadening of the spectrum when compared with the 
single-particle results. However, when calculating the ex¬ 
citonic binding energy Ex,b as the difference between the 
single-particle ground state transition energy and the ex- 
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citonic ground state transition energy, we find also large 
variations in Ex.t- The excitonic binding energy Ex,b 
is displayed in Fig. [5] as a function of the configura¬ 
tion number n. The calculated values scatter between 
36 meV and 90 meV. On average we find here an exci¬ 
tonic binding energy of approximately 56 meV (dashed 
line). The difference in the observed excitonic binding 
energies can be understood from the single-particle wave 
functions shown in Fig.[6l We find that for configuration 
20 the electron wave function, in comparison to configu¬ 
rations 5 and 75, shows very little charge density in the 
spatial region where the hole is localized, leading in the 
single-particle picture to a reduced wave function over¬ 
lap. From this one could expect a lower exciton binding 
energy in comparison to a configuration where electron 
and hole wave functions are localized in the same spatial 
region. As shown in Fig.[6l this is the situation for config¬ 
uration 75, which has consequently a very large excitonic 
binding energy (cf. Fig. [5]). The intermediate situation 
is realized for configuration 5 giving also an intermedi¬ 
ate excitonic binding energy. In general it is important 
to note that excitonic binding energies in nonpolar QWs 
are much larger than in c-plane systems. This effect is 
expected due to the absence of the macroscopic built-in 
field in nonpolar nitride-based QWsj ^^’^^ 

To shed even more light onto the impact of the 
Coulomb interaction on the wave functions, we use re¬ 
duced electron and hole density matrices to visualize the 
electron and hole densities under the influence of the 
Coulomb interaction. In general, the excitonic many- 
body wave function can be written as a linear com¬ 
bination of electron-hole basis states: 

Here |0) is the vacuum state, c^ the expansion coefficient, 
and ej (hj) denotes the electron (hole) creation operator. 
Electron and hole states are denoted by i and j, respec¬ 
tively. We can then define reduced density matrices for 
electrons and holes. For instance, for the electrons the 
density operator p® is given by: 

/S'" = X! = X! ■ (3) 

i,i' j i,i' 

The corresponding electron and hole densities are given 
by p® = (R|p®|R) and p^ = (R|p^|R), respectively. Fig¬ 
ure |9] depicts the calculated electron p® and hole p^ den¬ 
sities for the configurations 5, 20, and 75. We can infer 
from Fig. [S] that while the hole charge density is almost 
unchanged in comparison to the single-particle descrip¬ 
tion (cf. Fig. [6|), the electron charge density is signif¬ 
icantly affected by the attractive Coulomb interaction 
between electron and hole. The electron charge den¬ 
sity, under the influence of the Coulomb interaction, lo¬ 
calizes about the hole for all configurations. Thus our 
atomistic model predicts exciton localization in nonpo¬ 
lar (In,Ga)N/GaN QWs. This behavior is consistent 




FIG. 9: (Color online) Gronnd state electron (red) and hole 
(green) charge densities with Coulomb effects included for 
configurations 5 (Config 5), 20 (Config 20), and 75 (Config 
75). Results are shown for different view points. Light (dark) 
isosurfaces correspond to 5% (25%) of the maximum charge 
density value. Dashed lines indicate the QW interfaces. 
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FIG. 10: (Color online) Calculated excitonic ground state 
emission spectrum with light polarization vectors parallel (E || 
c) and perpendicular (E _L c) to the c axis. The (black) 
dashed-dotted line shows the experimental (unpolarized) PL 
spectrum. 


FIG. 11: (Color online) Calculated excitonic emission spec¬ 
trum for the light polarization perpendicular (red solid line) 
and parallel (blue dashed line) to the c axis. The spectral 
dependence of the corresponding DOLP, calculated according 
to Eq. 0, is shown by the (green) dashed-dotted line. 


with the experimental observation of single-exponential 
PL decay transients discussed in detail in Sec. IIVBI and 
observed experimentally by other groupsi ^^:^° Also the 
strong wave function overlap independent of the hole lo¬ 
calization is consistent with the measured decay times 
being constant over the spectrum. 

To further compare the theoretical results with the 
experimental PL data, we have calculated the excitonic 
ground state emission spectrum following Ref. [s^ Since 
the experiment indicates a large DOLP (> 90%) we 
have performed calculations for two different light po¬ 
larization vectors Gp^i, namely = (1,0,0)"^ and 

Gp ii = (0,0,1)"^. This means that the selected electric 
field E is perpendicular to the c axis (ep_^) and parallel to 
the c axis (e^ y), respectively. In the calculations we have 
assumed growth along the y axis; therefore the chosen 
light polarization vectors reflect the experimental setup 
described in Sec. IIV Al The calculated excitonic ground 
state emission spectrum is shown in Fig. [10] for (red 
solid line) and ep y (blue dashed line) together with the 
experimental (unpolarized) PL emission spectrum (black 
dashed-dotted line). Several different features are visible 
in Fig. Uni 

First of all we find a good agreement between the calcu¬ 
lated (solid red) and the experimentally determined emis¬ 
sion energy (black dashed-dotted line), given the slight 
uncertainties in In content and well width as discussed in 
Sec. lIII Al Second, we find also theoretically a very broad 
emission spectrum. For E _L c (op^x), the theoretically 
determined FWHM is approximately lOI meV. Experi¬ 
mentally we find a FWHM value of 135 meV. Different 
factors might contribute to the observed differences be¬ 
tween theory and experiment. For example, even though 
75 different microscopic structures may appear a large 


number, it could be the case that even more configura¬ 
tions have to be considered to fully resolve the measured 
FWHM, bearing in mind the large variations between 
different microscopic configurations (cf. Fig. |5|). Addi¬ 
tionally, if subtle nonrandom clustering effects exist, they 
may contribute to the broadening of the PL linewidth. 
However, the theoretically determined value of 101 meV 
for the FWHM is in reasonable agreement with the ex¬ 
perimental data. Third, Fig. (TOj shows that there is a 
large difference in the calculated intensities for ep_^ and 
Gp ||. In the theoretical analysis, the intensities are nor¬ 
malized to the intensity of Gp j_. If we calculate the DOLP 
from our theoretical spectrum via Eq. m based on the 
maximum intensities for Gp^^ and Gp y, we find a value of 
approximately 87%, which is slightly smaller than the ex¬ 
perimental values (> 90%). To compare our theoretical 
results in more detail with the experiment, we have also 
calculated the spectral dependence of the DOLP from 
our theoretical emission spectra for the two light polar¬ 
ization configurations depicted in Fig. (TOj Using Eq. o, 
our theoretical results for the spectral dependence of the 
DOLP are shown in Fig. [11] Similarly to Fig. ]31 the 
DOLP (green dashed-dotted line) is shown together with 
the excitonic emission spectrum for E _L c (red solid line) 
and E II c (blue dashed line). When comparing our theo¬ 
retical data with the experimental data shown in Fig. 0 
the theoretical results show slightly lower values than the 
experiment, plus the calculated DOLP is not as constant 
as the experimental data across the spectrum. Again, 
even though 75 configurations may appear a large num¬ 
ber, we show below that some of the structure in the 
DOLP spectrum of Fig. [TT]is due to a small number of 
exceptional states; more configurations would be required 
to reliably treat the importance of such states. This is 
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beyond the scope of the present study, since the present 
analysis gives already, in general, a good description of 
the experimentally observed spectral dependence of the 
DOLP. 

All in all the theoretical calculations indicate that ran¬ 
dom alloy fluctuations, explicitly taken into account in 
our model, affect the DOLP only slightly. The origin 
of the calculated high DOLP can be further understood 
by looking at the orbital character of the hole ground 
state/VBE state. The outcome of such an analysis is 
displayed in Fig. [TH where the orbital contributions to 
the VBE are shown as a function of the configuration 
number n. From a continuum-based calculation, neglect¬ 
ing the weak spin-orbit coupling, one would expect, due 
to the differences in the valence band effective masses 
along the growth direction and the positive crystal field 
splitting energy in GaN and InN, that the VBE is dom¬ 
inated by a single-orbital type {px- or py-like orbitals). 
Obviously such an analysis neglects the effects of alloy 
fluctuations. However, we can infer from Fig. [T^ that 
the VBE state in the different microscopic configurations 
is mainly made up of contributions from a single orbital 
type, in this case from px-V^e orbitals. Thus, for the or¬ 
bital character of the hole ground state the microscopic 
configuration is, in general, of secondary importance, al¬ 
though we note that there is an enhanced pz-like char¬ 
acter and a very low Pa;-like character in about 10%-20% 
of the structures studied. The dominance of the Px char¬ 
acter of the VBE explains the calculated high DOLP, 
which is in good agreement with the experimental data. 
Furthermore, as we know from Fig. [5] the hole ground 
state energies vary significantly between different config¬ 
urations. This gives rise to the broad emission spectrum 
shown in Fig. [TOl Since Fig. [T^] reveals that the orbital 
contribution to the VBE state is for the most part inde¬ 
pendent of the configuration number n, all these findings 
in combination explain why we observe only a weak spec¬ 
tral dependence of the DOLP displayed in Fig.[TTl in line 
with the experimental results displayed in Fig. [4] 


V. CONCLUSION 

In summary, we have presented a detailed experimental 
and theoretical analysis of the structural, electronic, and 
optical properties of m-plane (In,Ga)N/GaN QWs. The 
structural characterization using XRD, STEM, and APT 
provide us with consistent parameters for the QW com¬ 
position and dimensions and suggest that we may rea¬ 
sonably model the (In,Ga)N as a random alloy, although 
we are not able to rule out the possibility of some sub¬ 
tle deviations of the In distribution from randomness. In 
addition to the structural characterization, we used PL, 
P-PLE, and time-resolved PL measurements to analyze 
the optical properties of the system under consideration. 
Our results show a high DOLP, single-exponential PL de¬ 
cay transients, strong Stokes shifts, and a very broad PL 
linewidth. The strong Stokes shift and the broad PL line 



Config n 


FIG. 12: (Color online) Orbital contributions to the VBE as 
a function of the n different microscopic configurations con¬ 
sidered here. 


width are indicative of strong carrier localization effects, 
while single-exponential PL decay transients in nonpolar 
(In,Ga)N QWs are usually attributed to exciton localiza¬ 
tion effects. To shed further light on the experimentally 
observed properties, we have employed an atomistic TB 
model to achieve a microscopic description of the elec¬ 
tronic structure of the m-plane system including, on an 
atomistic level, strain and built-in field variations arising 
from the considered random alloy fluctuations. To be 
able to compare with the measured optical spectra, our 
model includes also excitonic effects via the GI scheme. 
The results of our calculations reveal strong hole wave 
function localization effects originating from random al¬ 
loy fluctuations. We find that alloy fluctuations have a 
much weaker effect on the electron wave functions. The 
observed localization effects lead to a signihcant broad¬ 
ening of the single-particle ground state energies and 
therefore the corresponding transition energies. When 
including Goulomb effects in the calculations we observe 
compared to c-plane systems strongly increased excitonic 
binding energies. This can be attributed to the absence 
of the macroscopic built-in field in the nonpolar m-plane 
system studied here. Coulomb effects lead also to the sit¬ 
uation where the electron wave function localizes about 
the hole. In other words, the theoretical calculations 
predict exciton localization effects, which are consistent 
with the single-exponential decay of the PL transients 
observed in the experiment. Moreover, when calculating 
the excitonic ground state emission spectrum we find a 
good agreement with the experimental PL spectra both 
in terms of FWHM and PL peak position energy. Addi¬ 
tionally, when calculating the DOLP, the theoretical data 
are in good agreement with the experimental result of a 
very high DOLP. Our findings indicate that the DOLP 
is, in general, only slightly affected by random alloy fluc¬ 
tuations. However, since the observed PL line width is 
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larger than the measured valence band splitting, wave 
function localization effects due to random alloy fluctua¬ 
tions are required to explain the spectral independence of 
the DOLP, given that the theoretical analysis shows that 
these localized states are dominated by a single orbital 
type. Overall, this combined experimental and theoreti¬ 
cal analysis provides clear insight into the basic physical 
properties of m-plane (In,Ga)N/GaN QWs. 
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